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ABSTRACT 

We construct mock galaxy catalogues to analyse clustering properties of a A cold 
dark matter (ACDM) universe within a cosmological dark matter simulation of suffi- 
cient resolution to resolve structure down to the scale of dwarfs. We show that there is 
a strong age-clustering correlation for objects likely to host luminous galaxies, which 
includes the satellite halo (subhalo) population. Older mock galaxies are significantly 
more clustered in our catalog, which consists of satellite haloes as well as the central 
peaks of discrete haloes, selected solely by peak circular velocity. This age dependence 
is caused mainly by the age-clustering relation for discrete haloes, recently found by 
Gao et al. , acting mostly on field members, combined with the tendency for older 
mock galaxies to lie within groups and clusters, where galaxy clustering is enhanced. 
Our results suggest that the clustering age dependence is manifested in real galaxies. 
At small scales (less than ~5 /i _1 Mpc), the very simple assumption that galaxy colour 
depends solely on halo age is inconsistent with the strength of the observed cluster- 
ing colour trends, where red galaxies become increasingly more clustered than blue 
galaxies toward smaller scales, suggesting that luminosity weighted galaxy ages do not 
closely trace the assembly epoch of their dark matter hosts. The age dependence is 
present but is much weaker for satellite haloes lying within groups and clusters than 
for the global population. 

Key words: galaxies: haloes - galaxies: formation - methods: N-body simulations - 
cosmology: theory - cosmology: dark matter 



1 INTRODUCTION 

A critical test of the ACDM model is whether it accurately 
predicts the clustering properties of galaxies formed within 
dark matter haloes and "subhaloes" that are satellite clumps 
withing larger host haloes. Subhaloes serve as hosts for vis- 
ible galaxies within clusters, groups, or larger galaxies, and 
so provide a natural basis for constructing simulated mock 
catalogs, whose clustering properties can then be compared 
with observed galaxies. 

Clustering of haloes depends on halo age, a phenomenon 
recently measured in ACDM simulations by Gao, Springel, & 
White (2005), who found that older haloes are more strongly 
clustered than younger haloes. The likely explanation is that 
haloes of a given mass generally form earlier within denser 
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regions. Thus, older haloes tend to populate denser regions, 
which naturally leads to stronger clustering with halo age. 
Sheth & Tormen (2004) measured such a trend between 
mean halo formation epoch and local over-density in simula- 
tion data, which was recently confirmed in larger simulations 
by Harker et al. (2005). The precise physical origin, however, 
of the age dependence is a subject of recent debate (see dis- 
cussion by e.g. Zentner 2006). Wang, Mo & Jing (2007), 
use numerical simulations to suggest that accretion onto low 
mass haloes in high density regions is inhibited by competi- 
tion with massive neighbours via tidal interactions and local 
dynamical heating, creating a correlation with halo age and 
environment. Sandvik et al. (2006) suggest that the forma- 
tion history and the current epoch environment of low mass 
haloes may be affected by their presence in massive pan- 
cakes and filaments at high redshift. The clustering age de- 
pendence, which has been confirmed by a number of authors 
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(Harker et al. 2005; Wechsler et al. 2006; Gao & White 2006; 
Jing, Suto & Mo 2006; Wang, Mo & Jing 2007), is strong 
for haloes of masses of 10 n ~ 12 /i _1 M® (Gao et al. 2005, 
Wechsler et al. 2006), which are likely to host galaxies, and 
decreases with halo mass, becoming insignificant for haloes 
more massive than ~10 13 /i _1 M (Gao et al. 2005). There 
is evidence that the age-clustering dependence may reverse 
sign for larger haloes (Wetzel et al. 2006; Jing, Suto & Mo 
2006; Gao & White 2006; Zentner 2006). 

Previous studies have focussed on the age dependence 
of clustering of discrete virialized haloes, and did not con- 
sider directly the contribution of satellite populations to the 
age dependence of clustering. Because a large fraction of 
galaxies belongs to groups and clusters, the clustering of 
the general galaxy population could have a strong depen- 
dence on subhalo ages. Galaxies of a given circular velocity 
will have formed earlier if they lie in groups or clusters to- 
day. Thus, we can expect that the contribution of group and 
cluster members will increase the tendency for older objects 
to be more strongly clustered. Also, recent studies relating 
subhalo numbers and distribution to age and to host halo 
properties hint that subhalo clustering could depend on sub- 
halo age. Recent simulations have shown that older haloes 
tend to host fewer subhaloes (e.g. Gao et al. 2004; Zentner et 
al. 2005; Taylor & Babul 2005; Zhu et al. 2006). Additionally, 
the clustering strength of virialized haloes is correlated with 
the numbers of their satellite haloes (Wechsler et al. 2006). 
Furthermore, subhaloes tend to lie nearer their host centres 
if they were either formed earlier (Willman et al. 2004) or 
were accreted earlier (Gao et al. 2004; Taylor & Babul 2005; 
see however Moore, Diemand, & Stadel 2004). 

In order to understand more fully the age dependence 
of subhalo clustering and its potential effects on observable 
galaxies, we analyze the relation between age and clustering 
within halo catalogs that include both the satellite haloes 
that populate group and cluster haloes as well as the dis- 
crete virialized haloes likely to host only a single galaxy. 
We construct a simple mock galaxy catalog wherein haloes 
are selected by peak circular velocity to roughly match the 
galaxy luminosities and abundances in large surveys. Our 
catalog is selected from a high resolution dark matter sim- 
ulation that resolves structures within a cosmological vol- 
ume down to the scale of dwarf galaxies. We stress that we 
are not attempting to create a realistic catalog of "simu- 
lated galaxies", but rather that we are merely using obser- 
vationally relevant circular velocities as a convenient means 
of assessing the potential dependence of clustering on age 
of haloes+subhaloes over a range that has the potential to 
host galaxies in the ACDM model. In § 2, we describe the 
simulations and the construction of the halo+subhalo cata- 
log. In § 3, we detail the age dependence of clustering in our 
mock catalog, the implications of which we discuss in § 4. 



2 NUMERICAL TECHNIQUES 
2.1 the simulations 

We use the parallel k-D (balanced binary) Tree (Bent- 
ley 1975) gravity solver PKDGRAV (Stadel 2001; Wadsley, 
Stadel & Quinn 2004) to model a 50 h~ Mpc cube, consist- 
ing of 432 3 dark matter particles of equal mass (the CUBEHI 



run of Reed et al. (2003; 2005ab). By modelling a relatively 
small cosmological volume, we are able to probe down to the 
small masses needed to resolve satellites within groups. The 
particle mass is 1.3x 10 s h~ 1 M ( 7 ) . A starting redshift of 69 
and a force softening of 5 /i _1 kpc (comoving) are used. This 
run adopts a ACDM cosmology with Q m — 0.3 and A 
f ) , 7 and the initial density power spectrum is normalised to 
a s = 1.0, consistent with WMAP (e.g. Bennett et al. 2003; 
Spergel et al. 2003). We use a Hubble constant of h = 0.7, in 
units of 100 km s _1 Mpc -1 , and assume no tilt (i.e. a primor- 
dial spectral index of 1). To set the initial conditions, we use 
the Bardeen et al. (1986) transfer function with T — fl m x h. 

2.2 mock catalog construction 

Mock galaxies are chosen from a catalog of haloes 
selected by circular velocity using the Spline Kernel In- 
terpolative DENMAX (SKID) halo finder (Stadel 2001; 
|http://wvw-hpcc7a stro. Washingto n. edu/tools/skid.htmI[ |. 
SKID haloes are identified using local density maxima 
to identify bound mass concentrations independently of 
environment. Note that SKID identifies discrete virialized 
haloes as well as subhaloes (satellite haloes). The radial 
extent of each SKID halo is determined by the distribution 
of bound particles, and no predetermined subhalo shape 
is imposed. The peak circular velocity of each subhalo, 
v c , P eak, is computed from the peak of the rotation curve 
v c (r) = (GM(< r)/r) 0,5 . The formation epoch is defined as 
the time at which v CtPea k of the largest progenitor (amongst 
all branches of the merger tree at a given epoch) reaches 
75% of its maximum value. A progenitor is defined as a 
halo with at least 30% of its particles incorporated into its 
descendent. Further detail on formation and accretion times 
of subhaloes can be found in a number of prior studies 
(e.g. De Lucia et al. 2004; Gao et al. 2004). 

The mock galaxies are selected to have a magnitude 
range similar to that of the Sloan Digital Sky Survey (SDSS) 
sample analyzed by Zehavi et al. (2002), -22 > M r > -19, 
though our results are not sensitive to the precise range. 
The absolute r-band magnitude M r of each SKID halo is 
estimated by applying v c ^ pea k to the Tully & Pierce (2000) 
variant of the Tully-Fisher (Tully & Fisher 1977) relation: 

M R = -21.12 - 7.65(logW R - 2.5), (1) 

where the linewidth W r is approximately twice v CyPea k 
(Tully & Fouque 1985). We select approximately 6,000 mock 
"galaxies" with 84 km s" 1 < v c ^ pea k < 206 km a -1 . In 
practice, our faintest simulated galaxies are SKID haloes of 
several hundred particles. While the Tully-Fisher magnitude 
assignment is subject to a number of uncertainties, including 
the fact that we apply this to ellipticals as well as spirals (see 
e.g. Desai et al. 2004), it provides a convenient method for 
building a catalog of mock galaxies with magnitudes compa- 
rable to those in galaxy surveys. The spatial abundance of 
the mock catalog is 4.8 x 10 -2 /i 3 Mpc -3 , which is 2.6 times 
that of the SDSS sample selected from the same magnitude 
range. However, we stress that our results are not sensitive 
to the abundance or to the limits used for inclusion into the 
mock catalog; i.e. , we are able to detect a clustering-age 
dependence for a range of u CiPea fc-selected catalogs in ad- 
dition to the one presented here, as we show later. Thus, 
even though the objects in our catalog are not expected to 
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describe precisely the galaxy population, we can still cap- 
ture many of the important clustering properties of the dark 
hosts of galaxies. 



2.3 correlation functions 

The spatial pairwise correlation function of galaxies is an im- 
portant cosmological test, as it quantitatively measures ba- 
sic clustering properties (see e.g. Peebles 1980). The spatial 
correlation function is calculated using the direct estimator 
(as in e.g. Governato et al. 1999): 



2N p (r) 
nlV(8V) 



(2) 



where N p (r) is the number of pairs in radial bins of volume 
SV, centred at r; n c is the mean space density of the catalog; 
and V is the volume of the simulation. We take into account 
our periodic boundary conditions when finding pairs. The 
correlation function is often approximated by a simple power 
law: 
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with £(ro) = 1, where ro is the correlation length. The rela- 
tive clustering amplitude between haloes and the mass dis- 
tribution is referred to as bias: 



(4) 



b 2 _ £,halo-halo{M , r, z) 

£d m {r,z)) 

For all error estimates, we use la poisson errors (equal to the 
square root of the number of pairs in each bin), which are 
likely to underestimate the true errors because they do not 
take into account clustering and sample variance (e.g. Croft 
& Efstathiou 1994). However, because we are interested pri- 
marily in the relative clustering between age-selected ob- 
jects, and not the true clustering strength, poisson errors 
are adequate for this study. 



3 RESULTS 

3.1 correlation functions of mass and haloes 

In Fig. [T] we plot £(r) and the bias factor (b(r) = 

\J ^haloes (r) /£,mass (r)) for SKID haloes with 
100, 150, and 200 km s _1 . Larger SKID haloes have steeper 
correlation functions and larger correlation lengths. The 
largest SKID haloes are "antibiased" (b < 1) with respect to 
the mass on small scales, but on large scales they are slightly 
more clustered than the mass. In general, haloes are "antib- 
iased" with respect to the mass, particularly on small scales 
for small haloes. This is consistent with small scale "antib- 
ias" found in previous simulations (e.g. Jenkins et al. 1998; 
Colin et al. 1999; Kravtsov & Klypin 1999; Yoshikawa et 
al. 2001; Diemand et al. 2004; Reed et al. 2005b) and may 
be caused by merging or destruction of subhaloes in high 
density regions (e.g. Jenkins et al. 1998; Klypin & Kravtsov 
1999). However, previous similar studies have found that 
the correlation function slope does not become shallower 
at small scales (Colin et al. 1999; Kravtsov et al. 2004; 
Neyrinck, Hamilton & Gnedin 2004; Conroy, Wechsler & 
Kravtsov 2006). The reasons for the difference are not clear, 
but we we note that the correlation function is a combination 



1000 



100 



^•10 




particles 
haloes > 200 km s~ 
haloes > 150 km s~ 



haloes > 100 km s 1 
haloes > 50 km s -1 

J I I I I '''I L_ 




r (h _1 Mpc) 



Figure 1. Top panel: The real space two-point correlation 
function for our SKID haloes with v c pea i c > 50, 100, 150, and 
200 km plotted along with §(r) for particles (mass). The 

horizontal line is £(r) = 1. Bottom panel: The halo bias b(r) = 



ihaloesir) / (.massir) . SKID haloes are selected by local density 
maxima, independent of environment, and include the centres of 
discrete virialized haloes in addition to self-bound satellite haloes 
(subhaloes). 



of a number of non-powerlaw components (central-satellite, 
central-central, and satellite-satellite) , so there is no a priori 
expectation that £(r) should follow a power-law (e.g. Benson 
et al. 2000; Berlind & Weinberg 2002; Kravtsov et al. 2004). 
In fact, the small scale departure from a power law that we 
see begins approximately where the satellite-satellite term 
begins to dominate £(r). At these scales, £(r) could be sen- 
sitive to a number of issues that affect the relative contribu- 
tions of these components. For example, the number of mas- 
sive clusters, which can dominate the satellite-satellite term, 
can be affected by run to run "sample variance" or box size 
(see e.g. Reed et al. 2007 and references therein). Differences 
in halo finder behavior at small scales could also be impor- 
tant. Further study is warranted, though our conclusions are 
not dependent on the smallest scales. It is difficult to quan- 
tify the precise scale below which the correlation function 
will no longer be robust. However, Reed et al. (2005b) in- 
dicate that the subhalo distribution is robust down to 100 
/i _1 kpc for simulations of similar resolution. For this reason, 
we have plotted all correlation functions only down to 100 
/i _1 kpc. 



3.2 the age dependence of the mock galaxy 
catalog correlation function 

We plot £ (r) for our mock galaxy sample in Fig. [2j binned 
according to formation times. Older catalog members are 
significantly more clustered for all pair separations. The old- 
est 10% is most preferentially clustered at small scales, with 
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Figure 2. The real space two-point correlation function £(r) for 
our mock catalog, which consists of v c iPea fc-selected SKID haloes 
sorted by age, determined by the time when the circular velocity, 
v c,peafei of the largest progenitor reaches 75% of its maximum 
value. The dot-dashed lines corresponds to £(r) = 1, and £(r) = 





r (h _1 Mpc) 

Figure 4. Same as Fig. [2] but for group members within the 
mock catalog, where a group is taken to be an friends-of-friends 
identified group of at least 3.2 X 10 13 Mq/i — Age percentiles are 
based on the group members only as opposed to from the full 
catalog. 



a clustering amplitude of ~10x that of the full mock cat- 
alog for separations less than ~1 h~ Mpc. The differences 
between the clustering of the young samples and the full 
catalog are smaller, but are significant. There is little dif- 
ference in the spatial correlations of the youngest 10% and 
the youngest 50%. The striking visual appearance of the 
age-clustering dependence is seen in Fig. [3] which shows the 
redshift zero simulation snapshot divided into the youngest 
and oldest 20% subsets of the mock catalog. 

Fig. [4] shows the correlation function for members of 
groups or clusters larger than 3.2 x 10 13 Mq/i _1 . There is 
some age dependence of the clustering of group members, 
but it is limited mainly to small pair separations, and is 
significantly weaker than that found for the full mock cata- 
log. Clustering of field objects, shown in Fig. [S] has a strong 
age dependence, though not as strong found in the complete 
sample. To determine group membership, groups are iden- 
tified using friends-of-friends (FOF; Press & Davis 1982; 
Davis et al. 1985), wherein the FOF haloes consist of par- 
ticles of separated by less than 0.2 times the mean inter- 
particle separation. The group extent is subsequently com- 
puted assuming a virial overdensity of approximately 100 
times the critical density (Eke, Cole, & Frenk 1996), and 
mock catalog members whose centre of mass lies within 
this region are assigned membership to that group. The 
overall correlation amplitude is a significantly higher for 
group members than for field members, an unsurprising re- 
sult given that group members are selected deliberately from 
within regions of high density, and belong to massive haloes, 
which are strongly clustered due to the well-known mass- 
clustering relation. 



The age-clustering relation in our mock catalog is likely 
due to a combination of causes. For the field sample, the ob- 
vious mechanism is the Gao et al. age-clustering correlation 
for discrete virialized haloes. For the full field plus group 
and cluster catalog (Fig. 0, the age dependence is stronger 
than that of the field sample alone (Fig. [5} because group 
and cluster members, which are found in highly clustered 
environments, tend to be old. Even though group and clus- 
ter members comprise only ~ 10% of the full sample, their 
contribution to the clustering age-dependence is significant 
due to the strong age correlation with environment. For ex- 
ample, 80% of our group and cluster members are older than 
the median mock galaxy age; and group and cluster mem- 
bers are 10 times more likely to belong to the 10% oldest 
subset than to the 10% youngest subset. 

On small scales (less than ~ l/i~ 1 Mpc), dynamical in- 
teractions become important for group and cluster members. 
Upon accretion onto a group or cluster halo, the subhalo will 
spiral in via dynamical friction, undergoing tidal stripping 
in the process. This leads to a subhalo distribution where 
centrally located group and cluster members were accreted 
earlier and are older. This is the likely cause of the small 
scale age-dependence within the group and cluster subsam- 
ple. Finally, for mock galaxies belonging to groups or cluster 
of similar mass, the age dependence of their host group- 
group clustering may produce some effect on the correlation 
function, but it should be mild because our group and clus- 
ter hosts are larger than the ~ 10 13 Mq/i -1 mass threshold 
above which the age-clustering relation of discrete haloes 
becomes weak (Gao et al. 2005). 
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(a) youngest 20% (b) full sample (c) oldest 20% 

Figure 3. Snapshot shows the mock catalog, composed of v c pe(1 fc-selected SKID haloes, divided into subsets by age, youngest 20% 
(left), full sample (centre), oldest 20% (right). Each snapshot shows the entire 50 h~ 1 Mpc volume. 
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Figure 5. Same as Fig. 12141 but for field (non-group) members 
within the mock catalog. 



3.3 can the age-clustering relation cause the 
observed colour-clustering dependence? 

In this section, we perform a simple test to determine 
whether the magnitude of the age-clustering dependence 
seen in simulations could be sufficient to account for the 
observed colour-clustering dependence under the simple as- 
sumption that halo age is a proxy for galaxy colour. Here, we 
have split the catalog into an "old" and a "young" subsam- 
ple with a 2:1 ratio of old to young haloes, which matches the 
ratio of red to blue galaxies of the Zehavi et al. (2002) SDSS 
sample. In Fig. ,[6] we plot the redshift space two-point corre- 
lation function for these catalog subsamples. The magnitude 
of the age dependence in the mock catalog sample is compa- 
rable to the clustering colour-dependence in SDSS for pair 



separations larger than ~5 ft _1 Mpc; see Fig. 11 of Zehavi 
et al. (2002) for the SDSS comparison. However, the age- 
effect is much weaker relative to the observed colour trends 
for smaller separations, and is insignificant (within the un- 
certainties) for pairs separated by 200 /i _1 kpc or less. The 
observed galaxy clustering colour-dependence, however, ex- 
tends down to ~ 100 ft _1 kpc (e.g. Zehavi et al. 2002; 2005, 
Madgwick et al. 2003; Li et al. 2006). The redshift space 
correlation function suggests that the age dependence has 
the potential to account for the observed colour-clustering 
trends only at large galaxy pair separations. This provides 
an independent argument that galaxy luminosity-weighted 
ages, indicated by colour, are different from the ages of the 
host dark matter subhaloes in which they lie. 

Qualitative differences between age selected simulated 
haloes and colour-selected observed haloes are also apparent 
in the projected-space two-point correlation function, 



2iVp(r) 
H 2 A(SA) 



(5) 



where w(r p ) is the pair excess over random with projected 
separation r p binned with area 8 A, and {i is the mean pro- 
jected density of haloes over the projected simulation area A 
(see e.g. Peebles 1980). In Fig.[7J the differences between the 
overall slope and small scale amplitude of w p (r p ) are rela- 
tively small between the young and old simulation samples. 
However, observed red galaxies have a much steeper and 
larger amplitude w p (r p ) than blue galaxies (see e.g. Zehavi 
et al. 2002, Fig. 13; Zehavi et al. 2005, Fig. 13) for a wide 
range of luminosity selected samples (e.g. Li et al. 2006). 

We note that our overall correlation amplitude is signif- 
icantly smaller than that of the SDSS sample. This is due, at 
least in part, to our finite box size, which means that large 
scale density fluctuations are not fully and accurately repre- 
sented (e.g. Bagla & Ray 2005; Sirko 2005; Power & Knebe 
2006; Reed et al. 2007), and should not be interpreted as an 
indication of a conflict with observations. A further contri- 
bution to our lower clustering amplitude may be the higher 
spatial abundance of our mock catalog, which implies that 
we are selecting smaller, and hence less strongly clustered, 
objects than in the SDSS sample. 
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Figure 6. The redshift space two-point correlation function £(r) 
for our mock galaxy catalog of v c pea j.-selected SKID haloes 
sorted into old and young subsamples based on the time at which 
the circular velocity of the largest progenitor reaches 75% of its 
maximum. 
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Figure 7. Same as Fig. \E\ but showing the projected two-point 
correlation function w p (r p ). 



A caveat here is that the mock catalog is inherently 
different from the SDSS sample. Because the correlation 
function of the full mock catalog is not a power law, as 
is generally observed in real galaxies, one should question 
whether the relative correlations of age-selected samples will 
display the same properties as real galaxies. A truly realistic 
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Figure 8. Same as Fig. [6] but showing w p (r p ) for a catalog se- 
lected by v C}Pea k > 150kms _1 . 



simulated galaxy sample would of course require modelling 
correctly all baryon physics, including star formation and 
feedback at small scales, a task that is not feasible at this 
time. However, it is prudent to consider how the details of 
our mock catalog construction could affect the measured 
age-clustering signal. To enable better statistics, our mock 
catalog was selected to contain smaller galaxies than the 
SDSS sample. It is clear from Fig. Q]that a sample of mock 
galaxies with higher v CtPea k results in a correlation function 
that is closer to a power law, better matching observations. 
We thus make a test to show whether the age-dependence 
of the correlation function is sensitive to v CyPea k (or abun- 
dance). In FigO we show that the clustering dependence on 
age has similar qualitative behavior, although the age de- 
pendence is somewhat weaker, for a sample of larger mock 
galaxies selected purely by v c , pe ak > 150 km s _1 . This indi- 
cates that the scale dependence of age-clustering relation in 
our simulation is relatively insensitive to our specific choice 



of v c , P eak range. 



3.4 age dependence of the pairwise velocity 
dispersion 

Galaxy peculiar velocities, as measured by the pairwise ve- 
locity dispersion, also provide a valuable probe of galaxy 
clustering, as well as providing an important component 
for dynamical probes of the dark energy equation of state 
(e.g. Governato et al. 1997; Baryshev, Chernin & Teerikorpi 
2001). The 1-D pairwise velocity dispersion, cr||(r), is the ve- 
locity dispersion for particle pairs in the direction parallel to 
the line of separation. In Fig. [9^,, we plot cry (r) for our mock 
catalog and for a random subsample of particles. Old mock 
galaxies have substantially higher pairwise velocities than 
young mock galaxies, as expected given their higher degree 
of spatial clustering. The old galaxy cii(r) is ~ 50 km s _1 
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Figure 9. The 1-D pairwise velocity dispersion along the line of 
separation, orii (r), is plotted for our i) CiPe[! fe-selected mock catalog, 
divided into a young and an old subset. 



"hotter" than the combined sample, and the young sample 
is "cooler" by up to 200 km s _1 . The lower pairwise veloc- 
ities for the mock catalog with respect to the dark matter 
particles is consistent with its spatial "antibias", shown in 

Fig. m 



3.5 age dependence in a mass-selected sample 

One important difference between our i! C;P <. a fc-serected haloes 
and the Gao et al. age dependence of clustering strength in 
friends-of-friends haloes is that our sample has a large range 
in masses whereas the Gao et al. study considered clustering 
at fixed mass. Because there is both a mass-age dependence 
and a mass-clustering dependence in CDM models, it is use- 
ful to consider what is the relation between age and cluster- 
ing for SKID haloes at fixed mass. It should be noted that 
there are nontrivial dependencies of SKID masses on envi- 
ronment; for example, less mass will be found to be gravita- 
tionally self-bound in high density environments due mainly 
to tidal stripping, but also affected to some degree by the 
addition of the external potential in the computation of self- 
bound mass. Thus, it is not obvious that there should be a 
similar age dependence among SKID haloes at fixed mass as 
there is for friends-of-friends haloes. We show in Fig.llOlthat 
the age dependence of clustering strength is indeed present 
for SKID haloes. This effect has a strong dependence on pair 
separation and on halo mass wherein the oldest 20% of the 
1.3 — 1.6 x 10 10 /i -1 Mo haloes have a clustering amplitude 
approximately 10 times larger than that of the the youngest 
20% at scales of 0.5 /i _1 Mpc. At 10 11 ft~ 1 M Q , the effect is 
much weaker, consistent with little or no age dependence, 
though our uncertainties are large in this higher mass range 
due to the smaller number of haloes. 
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Figure 10. The two-point correlation function £(r) of age- 
selected samples of SKID haloes, based on the time at which the 
largest progenitor reaches 50% of its peak mass. Solid lines show 
£(r) for all SKID haloes within each specified mass range. The 
horizontal and sloped dot-dashed lines correspond to §(r) = 1, 
and 



?( r ) = (§) i respectively. 



4 DISCUSSION 

We have shown that the clustering age dependence found 
for discrete friends-of-friends haloes by Gao et al. is also 
present in a mock galaxy catalog that consists of dark haloes 
plus satellites selected by v CtPea k to correspond to probable 
galaxy hosts. This provides strong evidence that galaxy clus- 
tering properties depend on the assembly epoch of the dark 
matter hosts for a population that includes field galaxies as 
well as galaxies within groups. 

It has been suggested by e.g. Gao et al. that the age de- 
pendence of the clustering amplitude may be a problem for 
the models of galaxy formation (e.g. Kauffmann, Nusser, & 
Steinmetz 1997; Benson et al. 2000; Wechsler et al. 2001) and 
the models of galaxy clustering in the halo model (e.g. Sel- 
jak 2000; Cooray & Sheth 2002; Berlind et al. 2003; van den 
Bosch, Yang, & Mo 2003) that assume that statistical galaxy 
properties depend only upon halo mass. The result that the 
age dependence on clustering is weaker among group mem- 
bers suggests that any impact on these models will likely be 
strongest for the low mass haloes that typically host 1 
luminous galaxy. 

Though the clustering age dependence is very strong for 
the extremes of the mock galaxy population, the age depen- 
dence is generally weaker, and has a different scale depen- 
dence than the colour dependence observed in recent large 
surveys. We expect that these qualitative differences are not 
sensitive to the precise selection criteria of the simulated 
sample. The different behavior of the simulated clustering 
age dependence suggests that luminosity-weighted galaxy 
ages, i.e. colour, do not trace halo age. An apparent ob- 
served lack of correlation of stellar ages and halo ages is 
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evident by the general trend that massive (luminous) galax- 
ies tend to have old (red) stellar populations, in apparent 
contradiction to the inverse mass-age relation present in hi- 
erarchical structure formation. This naively suggests that 
local environmental effects may have a strong influence on 
galaxy stellar populations. However, the different merger 
rates of haloes in regions of different density are also im- 
portant because progenitor haloes of a given mass will have 
been assembled earlier in more massive present day haloes, 
leading to earlier star formation (e.g. Mouri & Taniguchi 
2006; Neistein, van den Bosch & Dekel 2006). 

In any case, ages and colours of galaxies are expected 
to be influenced by a large number of astrophysical phenom- 
ena, including the suppression of star formation by winds, 
AGN, or ram pressure stripping (see e.g. Berlind et al. 2005; 
Bower et al. 2005). It may thus be difficult to decrypt a halo 
age-clustering signal in the real universe, as many of the 
influences on galaxy properties may depend on mass, envi- 
ronment, or other parameters (see e.g. Abbas & Sheth 2006; 
Cooray 2006). However, there may still exist some correla- 
tion between stellar population and halo age, if for example, 
major mergers trigger major starbursts, or if the age of the 
oldest stars in a galaxy is correlated with halo age. A re- 
cent semi-analytical study by Croton, Gao, & White (2007) 
suggests that that the clustering of group or cluster cen- 
tral galaxies should correlate with group host dark matter 
halo assembly age. Some evidence for an age-clustering trend 
of group haloes has recently been found by observing that 
groups of similar mass, whose central galaxies have more 
passive star formation, which may indicate earlier group as- 
sembly, are more strongly clustered (Yang, Mo & van den 
Bosch 2006). See however, an apparently opposite relation 
found by Berlind et al. (2006) , who find that massive groups 
tend to be less strongly clustered if they have redder cen- 
tral galaxies. To measure this effect in individual galaxies, 
one would need to measure accurately star formation his- 
tories of galaxies hosted by a narrow ranges of halo masses 
(see e.g. Heavens et al. 2004 for discussion of age measure- 
ments in SDSS stellar populations). Comparisons between 
the clustering age dependence among simulated haloes and 
the age dependence among galaxy stellar populations could 
then provide clues to the physics of galaxy formation. 

4.1 Summary 

Within a high resolution cosmological dark matter sim- 
ulation, we have examined the clustering properties of a 
mock galaxy catalog selected by v c , pea k to match approx- 
imately the luminosity range and number density of observ- 
able galaxies. 

• A strong clustering age dependence is found for mock 
galaxy catalogs that include both central haloes and satel- 
lite haloes, and is reflected in both spatial and in kinematic 
clustering measures. It is caused primarily by I) the age 
clustering relation for discrete virialized haloes, acting on 
field mock galaxies, and 2) the contribution of group and 
cluster members, which tend to be older, and are highly 
clustered due to their presence within massive dark matter 
hosts, thereby increasing the tendency for old members of 
the full sample to be highly clustered. 

• The strength of the clustering age dependence implies 
that it is likely to be manifested in real galaxies. The clus- 



tering age dependence is weaker than the clustering colour 
dependence in 2dF and SDSS for pair separations less than 
~5 /i _1 Mpc, and has a different scale dependence. This 
means that the observed clustering colour-dependence can- 
not be fully explained by assuming that stellar population 
ages trace halo ages. That is, one cannot simply assume that 
red galaxies lie in old haloes and blue galaxies lie in young 
haloes. The clustering colour dependence must be influenced 
by additional processes that affect the baryons. 

• The clustering age dependence is weaker among group 
and cluster mock catalog members than for the general 
galaxy population. 
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